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Abstract 

Comparing dilepton data from Au+Au to p+p collisions at RHIC energies PHENIX revealed 
two striking features: (i) a large excess of the dilepton yield at low mass and low pr in central 
Au+Au and (ii) a significant enhancement of direct real photon radiation. Both features can be 
interpreted as thermal radiation from the hot quark matter produced at RHIC. Unlike at lower 
beam energies, at RHIC the contribution from the hot hadronic phase alone seems insufficient to 
account for the data. This suggests large a contribution from a partonic phase. Both experimental 
and theoretical progress are required for more detailed conclusions. 



1. Introduction 

Electromagnetic radiation is regarded as an ideal tool to probe strongly interacting matter 
Experiments at SPS produced a wealth of data and inspired much theoretical work that advanced 
our understanding of meson properties in matter and the relation to chiral symmetry restoration 
10, Ht]. The PHENIX experiment continues the experimental investigation at RHIC in a new 
energy domain. In this review I report on the recent PHENIX results and discuss the comparison 
to model calculations. 

Figure [T] shows a schematic e^e" pair mass distribution as expected from p+p collisions at 
RHIC energies. The e^e^ pairs are either from pseudo-sealer or vector-meson decays, typically 
after the collision on the time scale of the strong interaction, or from hard-scattering processes 
like open heavy-flavor production, Drell-Yan pair production or quark-gluon (qg) Compton scat- 
tering, all occurring early in the collision. In addition to these sources quark matter produced in 
heavy ion collisions may emit thermal radiation, leave imprints from medium modifications of 
the meson properties, and modify the contribution from hard-scattering processes. 

Isolating the medium eff'ects from the expected sources is the experimental challenge and this 
challenge is two fold. First one needs to determine and subtract the pair background, which in 
central Au+Au collisions is more than a factor of hundred larger than the signal. This subtraction 
is very involved and dominates the systematic uncertainties. After many years of analysis effort 
this subtraction is well under control (see [4, 5J for details). The second challenge is to accurately 
determine the expected sources shown schematically in Fig. [T] PHENIX was able to determine 
these expected sources from data on meson production cross sections (see Fig.|2]and l'^) and on 
heavy-flavor production 12, M, measured in the same experiment. 

2. Reference data from p+p collisions 

Data from p+p collisions are used to test the analysis method and serves as precision refer- 
ence for data from heavy ion collisions. The data are published in [6]. They are compared to the 
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Figure 1: Schematic view of the e^e" pair mass distribu- 
tion from p+p collisions at 200 GeV. 
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Figure 2: Meson production cross-sections from 200 
GeV p+p collisions measured by PHENIX. The data for 
pseudo-sealer and vector mesons are compared to the dis- 
tributions derived by scaling the shape of the distributions 
by transverse mass and fitting the overall normalization to 
data d. 



expected yield in Fig. [3] It is evident from the ratio data-to-expected given in the lower panel 
that the data are in excellent agreement with the expectation within the systematic uncertainties 
of about 30%. 

The data are very precise and can be used to independently measure some of the expected 
contributions. This was done for the vector mesons u), (p [9] and II Oil as well as for heavy- 
flavor production |6] and qg Compton scattering. The later contribution can be isolated in the 
kinematic region where the mass is smaller than the transverse momentum ll Ull2ll . The ratio of 
direct-to-inclusive virtual photons can then be extrapolated to mass equal zero and compared to 
perturbative QCD calculations ifisll . As shown in Fig.|4]the agreement is very good. 



3. Data from Au+Au collisions 

With our understanding of the expected sources we have a solid base line to search for effects 
off the medium. Already from a first glance at the mass distribution from Au-hAu collisions 
(Fig E) one finds a striking enhancement of the data above the expected yield that was established 
in p-Hp collisions. In the mass range from 150 to 750 MeV/c^ the yield is a factor 3.4+0. 2(ifaf.) + 
1.3(i3'i.) ± Q).l{model) above the expectation. 

In order to investigate the enhancement in more detail one can study its centrality and mo- 
mentum dependence. The centrality dependence of the integrated yield is shown in Fig. |6] for 
two mass ranges m < 100 MeV/c^ and 150 to 750 MeV/c^, respectively. The yield in the lower 
mass bin is dominated by pairs from Dalitz-decays and consequently follows the centrality 
dependence of the pion yield. However, the yield in the higher mass range, which includes the 
enhancement, deviates from the expected dependence and indicates that central collisions con- 
tribute most of the observed enhancement. The yield from tttt or qq annihilation in the medium 




Figure 3: Electron-positron pair data from 200 GeV p+p collisions observed by the PHENIX experiment f^. The solid 
line depicts the contribution from hadron decays estimated on the basis of the most precise available data on neutral 
meson production and charm/bottom production. Data are connected for efficiency and require both electron and positron 
in the detector acceptance. 




Figure 4: Ratio of direct-to-inclusive dilepton Figure 5: Inclusive e+e" mass spectrum from Au+Au coUisions at ^/s 

pairs in the kinematic range m « pr- The data = 200 GeV measured by PHENIX. The presentation is identical to that 

are compared to perturbative QCD calculations in Fig. [3] The statistical errors are shown as bars and the systematic 

for direct photon production 1131 . en'ors are marked independently by a band around each data point. 
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Figure 6: Electron-positron pair yield per number of par- 
ticipants for two mass ranges sliown as function of colli- 
sion centrality. The data are compared to the expectation 
from known sources. 



Figure 7: Inverse slope of the nir distribution for mini- 
mum bias Au+Au collisions as a function of mass for two 
different regions for mr. Different from previous figures 
the data were corrected for pair acceptance. Again the 
data are compared to the expectation from known sources. 



should also increase faster than the number of participants, but it is not clear if the data can be 
explained quantitatively. 

Splitting the data in bins of pair pr reveals that in the mass range from 150 to 750 MeV/c^ 
the yield is enhanced at all pr, but not uniformly. The enhancement is largest at the lowest pr- 
To quantify this more PHENIX has measured the transverse mass (mj-) distributions, corrected 
for pair acceptance, and determined local slopes as function of pair mass. The inverse local slope 
or Tgff for the two regions < nir < 1 GeV and 1 < nir < 2 GeV is plotted in Fig. [T] as a 
function of mass. Also shown is Tgff expected from hadron decays alone. In general neither 
expectation from hadron decays nor the data show a strong variation of the effective temperature 
over the observed mass range. In the higher mj range the effective temperatures are similar to 
the expectation, 280 Me V compared to 220 MeV in the data. In the lower nir range the data have 
a significantly smaller slope of 120 MeV than the expected 200 MeV. Unfortunately, limited 
statistics prohibits a more detailed analysis. More precise data are needed to subtract the hadron 
decay contribution and isolate the enhancement. However, the effective temperature of the low 
niT component is a rather robust measurement and the result will not change much once the 
excess can be isolated. 

I now turn the attention to the contributions from hard scattering processes. Due to the 
limited statistical accuracy not much can be said about the contribution from charm production. 
The data shown in Fig. |5]are consistent with the extrapolation of the p-i-p cross-section using 
binary scaling, which is labeled cc (PYTHIA) in Fig. |5] However, data on charm production 
|@] indicate a large degree of interaction of charm quarks with the medium in central Au-i-Au 
collisions, which must also modify the angular correlation between c and c. If the momentum 
direction of the c and c quarks would randomized in the medium the data should be compared 
to the curve labeled cc (random correlation), which leaves room for other contributions above 1 
GeV/c2. 

Interesting results emerge when looking at the kinematic region of m « pr, where direct 
virtual photon emission was observed in p+p collisions. Details of the analysis, in particular 
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Figure 8: Ratio of direct-to-virtual photons. The plot is 
similar' to Fig.|4]but for Au+Au rather than p+p collisions. 
The pQCD prediction was scalled by the number of binary 
collisions. 




Figure 9: Invariant cross section of direct photon produc- 
tion deduced from Fig.[8] Results from p+p collisions are 
compared to pQCD calculations, while results for Au+Au 
are compared to a fit to p+p data scaled by the number of 
binary collisions. 



a discussion of the validity of extrapolating virtual photons back to the real photon point at 



trapc 

mass equal zero, can be found in II ll Il2ll . The final result for the direct-to-inclusive virtual 



photon ratio, shown in Fig. [8] indicates a significant excess beyond the expectation from initial 
state qg Compton scattering. Fig. |9] shows the direct photon spectrum resulting from Fig. [S] 
Here the virtual photon data were extrapolate to the real photon point, which in this case means 
setting the direct-to-inclusive ratio for virtual photons equal to the one for real photons. The data 
are consistent with results obtained from calorimeter measurements at higher pj in the same 
experiment. While p-i-p data follow the expectation from pQCD an excess is found in Au+Au. 
The inverse slope of the excess, after subtracting the scaled p-i-p data are Tgjf ~ 220 MeV. 
This effective temperature can be interpreted as Umit on the initial temperature of the reaction 
volume. Work to confirm these results with data from real photons is underway and has produced 
consistent upper limits. 



4. Model comparison and discussion 

Much work has gone into modeling dilepton and photon production from the fireball created 
in heavy ion collisions at SPS energies |2]. After some 20 years of experimental and theoreti- 
cal development the following picture has emerged. Electromagnetic radiation is predominantly 
emitted from hadron-hadron collisions at times when the spectral function of mesons are sig- 
nificantly broadened compared to vacuum. The colliding hadrons participate in the collective 
expansion and consequently the dileptons show radial flow. While emission at the partonic level 
is also observed it does not contribute significantly to the over all yield [2]. 

The same sources must also contribute to dilepton production at RHIC energies. In Fig.fTOl 



models, which describe SPS data, are confronted with dilepton mass distribution lll4ll5Lll6i1 . Ev- 

5 



idently, the predicted yield underestimates the data and these models are insufficient to describe 
RHIC data. Most likely the contributions from the partonic phase must be more significant than 
realized. At closer inspection it turns out that the model calculations shown in Fig. [10] include 
only ^^-annihilation but not the contribution from qg Compton scattering in the medium. 
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Figure 10: Comparison of Au+Au data to various models that successfully describe data from lower energies. Only 
calculation assuming broadening of the meson masses in the medium are shown. 
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Figure 11: Comparison of Au+Au data on direct photon emission with theoretical models. 
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Figure 12: Comparison of Au+Au data to various models that successfully describe data from lower energies. The 
labeling scheme is identical to Fig. 10. 



Model calculations of thermal photon emission based on perturbative QCD show that qg 
Compton scattering is very important 11211 . In a QGP the gluon densities might be so large that qg 
Compton scattering outshines ^^-annihilation and the perturbative approach may be insufficient 
all together. In Fig. [TT] various models are compared with real photon data [17]. All models agree 
roughly with the observed slope. The initial temperature in the models varies from 300 to 600 
MeV depending on the initial conditions, in particular the formation time. However, there are 
large differences in the predicted cross sections, which need to be addressed in future work. 

Last I return to dilepton spectra, but in a double differential representation in mass and nir- 
In Fig.[T2]the same calculations already shown in the mass projection (Fig.fTOli are compared to 
the rriT spectra. In the mass range from 500 to 750 MeV/c^ the models describe the data except 
for the lowest mj, as expected from Fig.[TO] At lower masses the models fall short at all mj- It 
is interesting to point out that the data above mr ~1 GeV/c^ were used to determine the direct 
photon yield. If qg Compton scattering would be consistently included in the models it would fill 
some of the difference at lower nir, though probably at the lowest mr it would not account for all 
data. Another important consequence may be that the models then would likely overestimate the 
yield in mass range 500 to 700 MeV/cl Obviously more complete theoretical calculations are 
required. In the highest shown mass range the models predict widely different yields, most likely 
due to different contributions from hard scattering processes. Again this needs to be looked into. 

To conclude this section I want to underline the importance of consistent model calculations 
that not only have a realistic description of the space time evolution and hadron spectra, but that 
also describe dilepton and photon emission simultaneously. 
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5. Summary and outlook 



Dilepton data analysis of the PHENIX experiment has matured significantly over the past 
years. Experimental issues, most importantly the subtraction of pair background, are under con- 
trol. Data from p+p collisions establish a precise reference to detect medium effects in Au+Au 
collisions. PHENIX has discovered such effects, namely a significant enhancement of dilepton 
production in mostly central Au+Au collisions. The enhancement covers the mass range from 
roughly 100 MeV to 750 MeV and is focused at low nir- Most strikingly it exhibits a promi- 
nent soft contribution with Tgff ~12Q MeV, which is independent of mass over the accessible 
range. PHENIX has also presented a first measurement of thermal photons indicating an initial 
temperature larger than 220 MeV. 

Though much work on theoretical models is needed to make final conclusions some firm 
conclusions can already be made. Thermal radiation from the hadronic phase, mostly from 
TTTT-annihilation with collision broadening has emerged as the main contribution to dilepton en- 
hancement observed at the SPS [2]. These contributions that obviously must exist also at RHIC 
energies remain insufficient to explain the PHENIX data. More complete calculations of dilepton 
radiation from the partonic phase need to be made. Model calculations of thermal real photon 
production is consistent with the data and hint towards initial temperatures larger than 300 MeV. 
However, big differences between various calculations remain. 

PHENIX has just started their dilepton program, more data from Cu+Cu , p+p , and d+Au 
collisions are on tape and are being analyzed. During RHIC Run 9 the collaboration successfully 
commissioned the hadron blind detector (HBD) and we can look forward to precision Au+Au 
data from PHENIX from Run 10. 
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